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INTRODUCTION 


The  presence  of  speckle  in  imagery  produced  with  coherent  illumi¬ 
nation  (e.g.,  synthetic  aperture  radar)  reduces  the  detectability  of 
objects  in  the  image%,^[-l-9] .  It  also  reduces  the  effectiveness  of 
some  computer  algorithms  (e.g.,  edge  detection)  designed  for  auto¬ 
matic  image  analysis. 

Thus  it  is  desirable  to  reduce  speckle  in  synthetic  aperture 
radar  (SAR)  images  both  to  assist  radar  image  interpreters  and  to 
preprocess  images  for  automatic  recognition  algorithms  on  computers. 
The  goal  is  to  smooth  out  the  speckle  but  at  the  same  time  preserve 
features  of  interest  like  edges,  strong  returns,  etc. 

vr 


SUMMARY 


2.1  THE  GEOMETRIC  FILTER 

A  non-linear  filter  was  developed  to  reduce  speckle  in  SAR 
imagery  while  at  the  same  time  preserving  spatial  information.  This 
filter  is  fundamentally  different  from  any  of  the  speckle  reduction 
techniques  that  have  appeared  in  the  literature  [10-33],  It  is  based 
on  geometric  concepts  ([34]  -  Appendix  A),  and  is  referred  to  as  tne 
geometric  filter. 

The  original  geometric  filter  was  designed  for  8-oit  imagery. 
It  is  an  iterative  algorithm  and  generally  about  four  to  ten  itera¬ 
tions  is  optimal  although  just  two  or  three  iterations  will  reduce 
the  speckle  significantly. 

This  8-bit  filter  was  originally  implemented  on  a  De  Anza  IP5500 
digital  video  processor  using  a  VAX  11/980  as  a  host  system.  Running 
time  was  about  14  s  per  iteration  for  512  x  512  images.  It  was  later 
implemented  on  a  4-stage  configuration  of  the  CYTO-HSS  High  Speed 
Cytocomputer  System  (recently  developed  at  ERIM  -  see  [46])  and  run¬ 
ning  time  was  reduced  to  0.2  s  per  iteration  for  images  of  the  same 
size.  This  hardware  can  also  be  configured  with  32  processor  stages 
which  would  further  reduce  the  execution  time  of  0.025  s  for  such 
images . 

The  geometric  filter  was  compared  favorably  ([34]  -  Appendix  A) 
with  the  3x3  median  filter  [21-29],  It  was  also  compared  favorably 
([35,  36]  -  Appendicies  B,  C)  with  look-averaging  [20,  32]. 

Exploration  of  filters  in  a  class  of  non-linear  filters  failed 
to  produce  any  that  performed  as  well  as  the  geometric  filter. 

2.2  THE  GENERALIZED  GEOMETRIC  FILTER 

One  iteration  of  the  geometric  filter  involves  32  basic  opera¬ 
tions,  16  of  which  increase  pixel  values  and  the  otner  16  decrease 
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pixel  values  (see  Appendix  A).  In  the  original  version  discussed 
above,  a  pixel  value  can  be  increased  or  decreased  by  at  most  one  by 
each  of  these  basic  operations.  This  original  version  was  developed 
for  use  on  8-bit  SAR  imagery. 

A  generalization  of  the  geometric  filter  was  developed  for  use 
on  16-bit  SAR  imagery.  Because  of  the  greater  dynamic  range,  many 
iterations  of  the  original  version  would  be  required  to  reduce  the 
speckle  significantly.  The  generalized  version  has  an  integer-valued 
control  parameter  n  which  allows  for  adjustment  to  the  dynamic  range 
of  the  image.  Each  of  the  32  basic  operations  referred  to  above  can 
increase  or  decrease  a  pixel  value  by  at  most  n.  With  n  =  1,  the 
operation  of  the  generalized  version  is  identical  to  that  of  the 
original  version.  A  program  that  computes  one  iteration  of  the  gen¬ 
eralized  geometric  filter  is  presented  in  Appendix  D. 

Experiments  were  performed  to  determine  optimal  settings  for  the 
control  parameter.  The  maximum  amount  a  pixel  value  can  be  changed 
in  a  single  iteration  is  16  n.  For  example,  if  there  is  a  pixel  in 

the  image  whose  value  is  more  than  16  n  greater  than  tne  value  of 

any  other  pixel  in  the  image  then  its  value  will  be  reduced  by  16  n. 
Thicker  spikes  will  be  reduced  more  slowly. 

The  speckle  index  (see  Appendix  A)  was  computed  before  and  after 
running  the  filter  with  various  settings  of  the  control  parameter. 
The  extent  to  which  the  spatial  information  was  preserved  was  judged 
subjectively. 

By  this  somewhat  ad  hoc  procedure  it  was  found  that  the  best 
results  were  achieved  by  halving  the  value  of  the  control  parameter 
on  each  successive  iteration.  It  was  also  found  that  the  maximum 
change  in  pixels  values  (=  16  n)  for  the  first  iteration  should  be 
about  one  quarter  of  the  maximum  pixel  value  in  the  image.  This 

leads  to  the  following  method  for  setting  the  control  parameter.  On 

the  first  iteration  the  control  parameter  is  set  equal  to  n^  wnere 
nl  is  the  smallest  power  of  two  which  is  greater  than  or  equal  to 
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maximum  pixel  value  in  the  image  divided  by  64.  The  control  param¬ 
eter  is  then  halved  on  each  successive  iteration.  On  the  last  itera¬ 
tion  it  has  the  value  1.  This  scheme  results  in  1  +  log^  n^ 

iterations.  If  a  16-bit  image  is  saturated  then  its  maximum  pixel 

value  is  2^  -  1  which  results  in  11  iterations. 

Table  2-1  shows  the  speckle  indices  and  maximum  pixel  values  for 
the  sequence  of  images  obtained  using  this  scheme  on  a  16-bit  512  x 
512  X-band  SAR  magnitude  image  of  an  area  in  North  Carolina  near 

Rayleigh-Ourham.  (This  imagery  was  made  under  USGS  Contract  No. 
14-08-0001-21748).  Table  2-1  also  shows  tne  reductions  in  tne  maxi¬ 
mum  pixel  value  for  each  iteration  and  the  maximum  possible  reauc¬ 
tions.  The  maximum  pixel  value  in  the  original  image  is  1946  and 
1946/64  =  30.4.  Therefore  the  control  parameter  settings  are  32, 
16,  8,  4,  2  and  1  for  the  six  iterations  performed.  (See  Figures 

2-1,  2-2  and  2-3.) 

2.3  CONTRAST  ENHANCEMENT  ALGORITHM 

Most  digital  display  devices  display  8-bit  images,  which  is  rea¬ 
sonable  since  the  human  vision  system  can  perceive  only  about  8  bits' 
worth  of  grey-levels.  In  order  to  display  a  16-bit  image  it  must 
first  be  scaled  down  to  8-bits.  The  simplest  method  is  to  use  a 
linear  scaling  wnich  maps  the  maximum  pixel  value  into  255  and  the 
minimum  into  0.  However,  if  the  image  contains  large  sections  with 
a  relatively  small  dynamic  range  plus  a  few  very  bright  points,  this 
method  will  result  in  using  very  few  quantization  levels  in  large 
sections  of  the  scaled-down  image  just  to  accommodate  the  few  bright 
points.  The  generalized  geometric  filter  produces  images  of  this 
sort.  8y  reducing  the  speckle  it  reduces  the  dynamic  range  in  large 
sections  of  the  image.  On  the  other  hand,  it  preserves  small  strong 
returns.  Although  the  pixel  values  of  these  returns  are  reduced, 
they  can  still  be  fairly  large. 


TABLE  2-1 

EFFECTS  OF  THE  GENERALIZED  GEOMETRIC  FILTER 


Iteration 

Con trol 
Parameter 
n 

Speckle 

Index 

Maximum 

Pixel 

Value 

Reduction 
in  Max 

Maximum 

Possible 

Reduction 

0 

N.A. 

.475 

1946 

0 

0 

1 

32 

.068 

1519 

427 

512 

2 

16 

.039 

1341 

178 

256 

3 

8 

.029 

1250 

91 

128 

4 

4 

.024 

1210 

40 

64 

5 

2 

.0215 

1190 

20 

32 

6 

1 

.0202 

1179 

11 

16 

Note!  the '"'<0th’*r- iteration  is  the  original  image. 
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The  contrast  enhancement  algorithm  is  basically  a  compromise  be¬ 
tween  the  global  linear  scaling  defined  above  and  an  adaptive  local 
linear  scaling  which  uses  local  maxima  and  minima.  If  the  adaptive 
local  linear  scaling  is  used  by  itself,  all  global  contrast  is  lost 
and  this  global  contrast  can  contain  useful  information.  For  exam¬ 
ple,  a  wooded  area  will  be  brighter  than  a  pasture  in  the  globally 
scaled  image. 


Let  f(m,n)  represent  the  filtered  image  and  let  G  be  its  global 
maximum  and  let  H  be  its  global  minimum.  Let  g^(m,n)  be  the  local 

maximum  of  f  over  a  N  x  N  square  centered  at  (m,n),  where  N  is  an 

odd  integer.  (If  the  pixel  (m,n)  is  close  enough  to  an  edge  so  that 
this  square  is  not  completely  contained  in  the  image,  the  local  max¬ 
imum  is  taken  over  the  part  of  the  square  which  is  in  the  image.) 
Let  h^(m,n)  be  similarly  defined  as  the  local  minimum  of  the  image 

f.  In  order  to  avoid  the  creation  of  artificial  edges  in  the  final 

image,  the  functions  g^  and  h^  are  smoothed  by  convolving  them 


with  the  function 


i,n)  =  | 


c  sinc  (inn:)"1  sinc(ir^T)n 


for  | ml  <  N 


0  otherwise. 


where  c  is  chosen  so  that 


ip(m,  n)  =  1, 


|  m  |  <  N 
n|  <  N 


We  define  g(m,n)  as  a  convex  combination  of  the  smoothed  local  maxi¬ 
mum  and  the  global  maximum  G; 


g(m,n)  =  8 


+  j,  n  +  k ) ip(  j  ,  k)  +  (1  -  8 )G,  (3) 


Ml  <  N 

Ik  1  <  N 


<J 
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where  s  is  an  input  parameter,  0  <  e  <  1.  (Appropriate  modifications 
are  made  near  edges.)  Another  function  h(m,  n)  is  defined  similarly 
as  a  convex  combination  (same  b)  of  the  smoothed  local  minimum  and 
the  global  minimum  H. 

Finally,  the  contrast  enhanced  image  is  defined  as 


f-j  (ro,  n) 


gfc-nj-  -~K(~m)~  n)  *  255  if  9("«.  ")  t  h(m,  n)  , 
~h~~ ~  ~  255  otherwise. 


The  image  in  Figure  2-1  was  obtained  by  simple  global  scaling 
(b  =  0)  of  the  original  16-bit  image  before  running  the  generalized 
geometric  filter  on  it.  A  16-bit  filter  image  was  produced  by  run¬ 
ning  six  iterations  of  the  filter  with  n  =  32,  16,  8,  4,  2  and  1, 
respectively.  Figure  2-2  is  the  8-bit  image  obtained  by  simple 
global  linear  scaling  (b  =  0)  of  the  filtered  image.  Figure  2-3  is 
the  8-bit  image  obtained  by  running  the  contrast  enhancement  algo¬ 
rithm  on  the  filtered  image  with  N  *  3  and  b  *  .95. 
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Geometric  filter  for  speckle  reduction 

Thomas  R.  Crimmins 


An  algorithm  is  described  which  reduces  speckle  noise  in  images.  It  is  a  nonlinear  algorithm  based  on  geo¬ 
metric  concepts.  Tests  were  performed  on  synthetic  aperture  radar  images  which  show  that  it  compares  fa¬ 
vorably  with  a  3  x  3  median  filter. 


I.  Introduction  and  Results 

It  is  desirable  to  reduce  speckle  noise  in  synthetic 
aperture  radar  images  both  to  assist  radar  image  in¬ 
terpreters  and  to  preprocess  images  for  automatic  rec¬ 
ognition  algorithms  on  computers.  The  goal  is  to 
smooth  the  speckle  but  at  the  same  time  preserve  fea¬ 
tures  of  interest  such  as  edges,  strong  returns,  etc. 

A  nonlinear  algorithm  based  on  geometric  concepts 
was  developed  to  accomplish  this.  It  is  an  iterative  al¬ 
gorithm  and  usually  about  ten  iterations  seems  to  be 
optimal  although  just  two  or  three  iterations  will  reduce 
the  speckle  significantly  (see  Table  I).  It  was  pro¬ 
grammed  on  a  DeAnza  IP5000  digital  video  processor 
using  a  VAX  11/780  as  a  host  system.  Running  time 
was  14  sec/iteration. 

Figure  1  is  a  synthetic  aperture  radar  image  of  Willow 
Run  Airport  in  southeastern  Michigan.  It  is  an  X -band 
strip  map  image  which  was  digitally  processed  on  the 
ERIM  digital  processor.  The  polarization  of  both 
transmitter  and  receiver  was  horizontal.  The  resolu¬ 
tion.  measured  as  the  half-power  width  of  the  impulse 
response  function,  is  6  m.  Figure  '2(a)  is  a  subimage  of 
Fig.  1.  The  pixel  spacing  in  Fig.  2(a)  is  3  m  and  the 
image  is  512  X  512  pixels. 

Figure  2(b)  shows  the  result  of  four  iterations  of  the 
geometric  Filter,  and  Fig.  2(c)  shows  the  result  of  ten  it¬ 
erations.  For  comparison.  Fig.  2(d)  shows  the  result  of 
applying  the  3x3  median  filter1-8  until  this  image  was 
obtained  which  is  of  period  2  under  the  median  F.lter. 
That  is.  one  application  of  the  median  filter  to  this 
image  will  change  it  (very  slightly)  but  it  is  invariant 
under  two  applications  of  the  median  filter.  (It  is  in- 
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teresting  to  note  that  this  means  that  this  image,  and 
hence  the  original  image,  have  no  median  roots.  This 
contrasts  with  the  1-D  case  in  which  median  roots  al¬ 
ways  exist.6)  The  3X3  median  filter  replaces  each 
pixel  value  with  the  median  of  the  nine  pixel  values  in 
its  3  X  3  neighborhood  window.  Figures  3(a).  (b).  and 
(c)  are  shaded  and  shadowed  versions  of  Figs.  2(a).  (c). 
and  (d),  respectively,  considered  as  2-D  surfaces  in  3-D 
space. 

Because  speckle’noise  is  multiplicative,  the  ratio  of 
its  deviation  to  its  mean  seems  to  be  a  reasonable  mea¬ 
sure  of  the  amount  of  speckle  noise  present  (see  Ref.  9. 
p.  25). 

Let  p(m,n),l  <  m,n  <  N,  be  the  pixel  values.  For  2 
<  m,n  <  N  -  1.  the  local  deviation  is  defined  as 

<r(m,n)  »  max  plm  +  i.n  +  />  -  min  plm  +  i.n  +  ;  >  ill 
- 1  <  i  <  1  - 1  <  I  <  1 

-1  <  i  <  1  -1 <  ;  <  1 

(This,  of  course,  is  not  the  true  deviation,  but  it  is  easy 
to  compute  and  we  believe  it  should  be  good  enough  for 
our  present  purpose.)  The  local  mean  is  defined  as 
l  i 

plm.nl*-  yj  ptm  +  i.n+;i  i2i 

9 1  j  ■  - 1 

Finally,  the  speckle  index  is  defined  as 

1  vj;i  cum.n) 

— - -  V - -  .  i3) 

(.V  -  2)2  plm.n  I 

Table  I  gives  the  speckle  index  for  the  original  image 
for  various  numbers  of  iterations  of  the  geometric  filter 
and  for  the  median  root.  Of  course  the  goal  here  is  to 
reduce  the  speckle  index  as  much  as  possible  while  still 
preserving  the  spatial  information  in  the  image.  The 
images  in  Figs.  2  and  3  indicate  how  well  the  spatial 
information  was  preserved.  Thus  the  speckle  index 
must  be  considered  together  with  the  corresponding 
images  in  order  to  determine  how  well  a  given  filter 
performs. 

The  geometric  filter  utilizes  a  technique  which  has 
a  history  starting  with  a  problem  which  has  nothing  to 
do  with  speckle  reduction.  The  original  problem  was 


Original  image 

Four  iterations  ot  geometric  filter 
Ten  iterations  of  geometric  filter 
Period  2  median  root 


mz 1 


Fig.  2.  lai  Original  image  ibi  Four  iterations  ot'  geometric  niter 
i c >  Ten  iterations  ot  geometric  filter  idi  Period  J  median  root  ■  > 
original  image. 


Fig  1.  Radar  image  of  Willow  Run  Airport.  The  boxed  section 
appears  in  Fig.  dial. 

to  generate  approximations  of  the  convex  hulls  of 
maximal  connected  subsets  of  the  foreground  of  a  bi¬ 
nary  image.  I A  binary  image  is  an  image  whose  pixels 
have  only  the  values  0  and  1.  The  foreground  is  the  set 
of  pixels  with  value  1.1  The  motivation  for  this  was 
enhancement  of  medical  imagery. 

II.  The  8-Hull  Algorithm 

The  convex  hull  of  a  set  is  the  intersection  of  all  the 
half-planes  containing  it.  An  approximation  to  this, 
called  the  8-hull,  is  defined  as  the  intersection  of  only 
those  half-planes  which  contain  the  set  and  whose  edges 
are  either  horizontal  or  vertical  or  lie  in  either  of  the  45° 
diagonal  directions.  The  8-hull  of  a  set  has.  at  most, 
eight  sides. 

F.  Lambeck  wrote  an  iterative  algorithm  (unpub¬ 
lished  i  which  generates  the  8-hull  of  a  set.  It  works  as 
ful lows  At  each  step  of  the  iteration,  the  value  of  a 
pixel  is  changed  from  a  0  to  a  1  if  its  neighboring  pixels 
have  ones  arranged  in  any  one  of  the  configurations 
<hown  in  Fig.  4.  The  blank  squares  can  be  either  zeros 
nr  ones.  If  enough  iterations  of  this  step  are  performed, 
eventually  the  8-hull  of  the  given  set  will  be  generated 
and  it  will  be  invariant  under  further  iterations.  An 
example  of  the  8- hull  of  a  set  being  formed  is  shown  in 
Figs.  o.  The  black  squares  represent  ones  and  the  white 
-quares  represent  zeros. 
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Fig.  3.  ia)  Shaded  and  shadowed  version  ot  original  image  'hi 
Shaded  and  shadowed  version  ot  ten  iterations  ot  geomeiru  filter  «  • 
Shaded  and  shadowed  version  of*  period  1  median  r<>. «t  <>l  original 


III.  Complementary  Hulling  Algorithm 

Next.  ■].  Gleason  developed  an  iterative  algorithm 
(unpublished)  for  smoothing  the  ragged  edges  ot  hinarx 
images  of  vehicles  obtained  by  slicing  gray-level  r.niar 
images.  We  will  refer  to  this  algorithm  as  the  comple¬ 
mentary  hulling  algorithm.  One  step  of  the  iteration 


15  May  1985  >  Vol  04.  No  10  APPLIED  OPTICS 


•409 


Fig.  6.  Allowable  90°  turns  in  the  boundary.  These  are  exceptions 
to  the  rule. 


(c)  <d) 

Fig.  5.  tat  Original  set.  t b)  One  iteration  of  8-hull  algorithm,  (c) 
Two  iterations,  (d)  Three  iterations.  This  is  invariant  under  further 
iterations. 

consists  of  the  following.  First,  one  step  of  the  8-hull 
algorithm  described  above  is  applied  to  the  set  and  then 
one  step  of  the  8-hull  algorithm  is  applied  to  its  com¬ 
plement.  In  other  words,  one  step  of  the  8-hull  algo¬ 
rithm  is  applied,  then  zeros  and  ones  are  interchanged, 
then  another  step  of  the  8-huli  algorithm  is  applied,  and 
finally,  zeros  and  ones  are  interchanged  again.  This  has 
the  effect  of  gradually  reducing  the  maximum  curvature 
of  the  boundary  of  the  set.  (Curvature  is  the  inverse  of 
the  radius  of  curvature.)  More  precisely,  with  a  few 
exceptions,  the  boundary  of  a  set  invariant  under  this 
algorithm  can  turn  a  maximum  of  45°  at  any  vertex. 
The  only  known  exceptions  to  this  rule  are  the  boundary 
segments  shown  in  Fig.  6  and  their  90°  rotations. 

Figure  7  shows  a  set  at  various  stages  of  the  iteration 
process.  The  ninth  iteration  is  invariant  under  further 
iterations.  Figure  8  shows  the  effect  of  this  algorithm 
on  a  set  with  a  tower  three  pixels  wide  protruding  from 
the  top.  If  one  imagines  going  up  the  left  side  of  this 
tower  and  then  down  the  right  side,  a  sharp  U-turn  must 
be  made  at  the  top.  This  U-turn  has  too  high  a  curva¬ 
ture  to  be  permissible.  Hence,  the  algorithm  reduces 
it  until  the  curvature  is  sufficiently  small.  A  tower 
seven  pixels  wide  is  shown  in  Fig.  9.  Here  the  U-turn 


Fig.  7.  (a l  Original  set.  (hi  One  iteration  ol  eomplement.irt  bulling 
algorithm  let  Two  iterations,  tdt  Three  iterations  <ei  Nine  iter, 
ations.  This  is  invariant  under  further  iteration' 
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Fig.  $.  la)  Original  set.  lb)  One  iteration  of  complementary  hulling 
algorithm.  ( c >  Two  iterations,  (d)  Fourteen  iterations,  (et  Fifteen 
iterations.  This  is  invariant  under  further  iterations. 


(c)  (d) 

Fig  '<  m>  Original  <et  i  hi  One  iteration  <it  complementary  hulling 
algorithm  u  <  I  wo  iterations.  <di  Three  iteration-  This  i-  in¬ 
variant  under  further  iterations. 


ran  be  made  without  requiring  too  high  a  curvature. 
Thu-,  the  turners  at  the  top  of  the  tower  are  rounded, 
but  the  tower  is  not  reduced.  It  is  this  reduction  of 
narrow  tower-  and  preservation  of  wide  towers  that 
-ugge-t  the  applu  abilit v  of  the  concepts  in  this  algo¬ 
rithm  to  -pet  kle  reduct  urn. 


IV.  Geometric  Filter 

A.  How  It  Works 

The  geometric  filter  uses  this  complementary  hulling 
technique  in  the  following  way.  First,  picture  the 
grav-level  (pixel  values  between  0  and  255  inclusive) 
radar  image  lying  horizontally  with  its  sides  running  in 
the  north-south  (N-S)  direction  and  its  top  and  bottom 
running  in  the  east-west  (E-W)  direction.  We  now 
construct  a  gray-level  surface  above  the  image  such  that 
its  height  above  any  pixel  is  proportional  to  the  value 
of  that  pixel.  Next,  this  surface  is  sliced  by  all  vertical 
planes  which  contain  a  line  of  pixels  running  in  the  E-W 
direction.  A  discrete  grid  is  defined  on  each  of  these 
vertical  planes  composed  of  vertical  lines  which  pass 
through  pixels  in  the  image  and  horizontal  lines  at 
heights  representing  the  integers  between  0  and  255 
inclusive.  The  points  in  this  vertical  grid  will  be  re¬ 
ferred  to  as  vertical  pixels.  Then  the  intersection  of  any 
one  of  the  vertical  planes  with  the  gray-level  surface  is 
a  curve  composed  of  vertical  pixels  in  the  vertical  grid 
[Fig.  10(a)).  Now  consider  the  umbra  of  this  curve. 
The  umbra  consists  of  all  vertical  pixels  in  the  vertical 
grid  on  or  below  the  curve  [Fig.  10(b)).  If  we  assign  the 
value  1  to  all  vertical  pixels  in  the  umbra  and  the  value 
0  to  all  other  vertical  pixels  in  the  vertical  grid,  the 
vertical  grid  can  be  considered  as  a  binary  image. 
(Nc<e:  It  is  not  necessary  for  the  algorithm  to  actually 
construct  this  binary  image.  We  construct  it  here  only 
for  pedagogic  purposes.)  Now  one  iterative  step  of  the 
complementary  hulling  algorithm  is  applied  to  the 
umbra,  which  is  the  foreground  of  the  vertical  binary 
image  just  constructed.  Recall  that  the  first  half  of  an 
iterative  step  of  the  complementary  hulling  algorithm 
consists  of  one  step  of  the  8-hull  algorithm.  Since  we 
are  only  concerned  with  the  line  forming  the  top 
boundary  of  the  umbra,  only  four  of  the  eight  configu¬ 
rations  shown  in  Fig.  4  need  he  used.  These  four  are 
shown  in  Fig.  4(a).  Similarly,  in  the  second  half  of  an 
iterative  step  of  the  complementary  hulling  algorithm, 
when  one  iterative  step  of  the  8-hull  algorithm  is  applied 
to  the  complement  [Fig.  10(c)],  only  the  line  forming  the 
bottom  boundary  of  the  complement  is  of  interest  and 
hence  only  the  four  configurations  shown  in  Fig.  4(b)  are 
used.  As  before,  the  complementary  hulling  algorithm 
step  is  completed  by  complementing  hack.  This  pro¬ 
cedure  is  performed  on  all  E-W  vertical  grids  simulta¬ 
neously  and  the  resulting  gray-level  image  replaces  the 
original  grav-level  image.  Now  the  same  procedure  is 
repeated  in  the  SW-NE  direction,  then  in  the  S  N  di¬ 
rection.  and  finally  in  the  SE-NW  direction.  This 
completes  one  iterative  step  of  the  geometric  filter. 

One  change  has  been  made  in  the  algorithm  described 
above  to  speed  the  speckle  reduction.  It  has  to  do  with 
the  application  of  an  iterative  step  of  the  .''-hull  algo¬ 
rithm  to  the  umbra  of  a  vertical  slice  as  described  above. 
Instead  of  changing  a  0  to  a  1  if  any  of  the  four  neigh¬ 
borhood  configurations  in  Fig.  4(a)  is  present,  these  tour 
configurations  are  used  -eparately  and  consecut ivel\ . 
That  is,  first  a  0  is  changed  to  a  1  ’f  the  first  ol  t  he-e  tour 
configurations  is  present.  Then  the  re-ult mg  image 
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Fig.  10.  (a)  Curve,  (b)  Umbra  of  curve,  (c)  Complement  of 

umbra. 


replaces  the  original  image.  Then  the  second  configu¬ 
ration  is  used.  etc.  Similarly,  the  four  configurations 
in  Fig.  4(b)  are  applied  separately  and  consecutively  to 
the  complement.  This  results  in  a  greater  amount  of 
modification  of  the  image  in  each  iterative  step  and 
hence  fewer  iterations  are  required.  It  also  causes  a 
greater  difference  between  the  reduction  rates  for  wide 
and  narrow  features  (see  next  subsection). 

B  Why  It  Works 

Speckle  has  a  wormy  appearance.  Considering  the 
image  as  a  curved  2-D  surface  in  3-D  space  ( as  described 
above),  speckle  appears  as  narrow  winding  walls  and 
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Table  II.  Reduction  Ratoa  lot  Square  Towor* 


Tower  sue 

Reduction  of 
maximum/iteration; 
number  of  vertical  pixel*, 
i.e..  unita  of  intensity 

1X1 

16 

2X2 

9 

3X3 

6 

4X4 

5 

5X5 

4 

6X6 

3 

7X7 

2 

8X8 

2 

Table  1*. 

Reduction  Rate*  tor  Wall* 

Reduction  of 

maximum/iteration; 

number  of  vertical  pixels. 

Wall  width 

i.e.,  units  of  intensity 

1 

12 

2 

6 

3 

3 

4 

3 

5 

2 

6 

1.5 

0 

valleys.  A  strong  target  return  appears  as  a  wide  tall 
tower.  The  entire  body  of  the  target  appears  as  an  even 
wider  high  plateau,  although  it  is  not  as  high  as  the 
tower  representing  a  strong  target  return.  The  geo¬ 
metric  filter,  through  iterative  repetition,  gradually 
tears  down  the  narrow  walls  and  fills  up  the  narrow- 
valleys.  It  also  tears  down  towers  and  high  plateaus, 
which  we  want  to  preserve.  However,  it  reduces  narrow- 
walls  and  valleys  faster  than  it  reduces  wide  towers,  and 
the  even  wider  plateaus  are  reduced  even  more  slowly 
than  the  wide  towers.  In  general,  the  wider  any  land¬ 
scape  feature  (valleys,  holes,  depressions,  walls,  towers, 
plateaus)  is,  the  more  slowly  it  is  reduced.  Thus,  only 
a  few  iterations  are  required  to  reduce  the  narrow 
speckle  walls  and  valleys  and  these  few  iterations  have 
very  little  effect  on  wide  high  towers  (hence,  strong 
target  returns  are  preserved)  or  on  the  even  wider  pla¬ 
teaus  (hence,  target  shape  is  preserved). 

Since  the  algorithm  works  on  both  the  umbra  of  the 
image  surface  and  its  complement,  it  is  essentially 
symmetric  with  respect  to  up  and  down.  Thus,  filling 
up  valleys  is  essentially  equivalent  to  tearing  down 
walls.  Tables  II  and  III  show  the  rates  at  which  the 
height  (maximum)  of  walls  and  square  towers  is  being 
reduced.  In  the  case  of  the  wider  towers  and  walls,  the 
reduction  rate  may  be  less  for  the  first  few  iterations. 
The  tables  show  the  rates  established  after  these  first 
few  iterations.  However,  this  lag  becomes  important 
for  wider  towers.  For  example,  the  maximum  of  a  14 
X  14  tower  is  not  reduced  at  all  until  the  fifth  itera¬ 
tion. 
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One  more  thing  should  be  mentioned  here.  Some¬ 
times  the  tower  representing  a  strong  target  return  is 
almost  as  narrow  as  the  speckle  wails  and  valleys. 
However,  such  a  tower  is  usually  much  higher  than  the 
speckle  walls  are  high  or  valleys  deep.  Thus,  although 
it  may  be  reduced  at  as  fast  a  rate  as  the  speckle  walls 
and  valleys,  the  walls  and  valleys  will  still  be  reduced 
before  the  tower  is  completely  torn  down  and  hence  the 
strong  target  return  will  still  be  preserved,  although  it 
will  be  dimmer  than  in  the  original  image. 

This  work  supported  in  part  by  the  Army  Research 
Office,  Physics  division,  under  contract  DAAG29-84- 
K-0204. 
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Geometric  filter  for  reducing  speckle 
Thomas  R.  Crimmins 

Radar  Division,  Environmental  Research  Institute  of  Michigan 
P.O.  Box  8618,  Ann  Arbor,  MI  48107 


Abstract 


A  run-linear  speckle  filter  based  on  geometric  concepts  is  defined  and  an  “xampl°  of 
its  e t i ec t 1 veness  on  synthetic  aperture  radar  imagery  is  shown.  A  comparison  with  look¬ 
averaging  is  made  using  artificial  imagery  with  synthetic  speckle. 

Introduction 


Tne  presence  of  speckle  in  imagery  produced  with  coherent  illumination  reduces  -he  !e- 
teccaoility  of  oojects  in  the  images^-  “  ®.  It  also  reduces  the  effectiveness  of  some 
computer  algortth  (e.g.,  edge  detection)  designed  for  automatic  image  analysis. 

Trie  geometric  filter  was  designed  to  reduce  speckle  in  synthetic  aperture  radar  (SAR) 
imagery  wnile  preserving  the  spatial  information  in  the  image  such  as  edges,  strnoo  re¬ 
turns,  etc.  it  is  an  iterative  non-linear  algorithm.  Usually  from  4  to  10  iterations  ire 
used.  '."ho  algorithm  was  implemented  on  a  DeAnza  IP  5500  digital  video  processor  using  a 
VAX  11  780  as  a  host  system.  Running  time  is  about  14  seconds  per  iteration  for  510  x  -10 
pixel  images. 


Definition  of  the  geometric  filter 

The  geometric  filter  algorithm  is  based  on  applying  a  single  iteration  of  an  iterative 
convex  .lulling  algorithm  alternately  to  the  image  and  to  its  complement  (negative  of  tne 
image).  For  details  of  the  geometric  derivation  and  an  intuitive  explanation  of  the  pro¬ 
gram  oeltw,  see  Crimmins9. 


IN 
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It  is  essentially  a  one-dimensional  algorithm  which  is  applied  successively  in  four  dif¬ 
ferent  urections  in  the  two-dimensional  image;  horizontal,  vertical,  and  the  two  diagonal 
directions.  Tne  algorithm  is  defined  as  follows.  The  image  is  a  function  ?(m,  D  where 
i  £  m  <  M,  i  S  n  S  N  and  the  values  of  f  are  integers  between  0  and  255.  A  border  of 
zeros  is  unded  to  the  image  so  that  f(m,  n)  is  defined  for  0  <  m  <  M  *  1  and  7  <  n  <  N  -  1 . 
An  auxiliary  image  g(m,  n)  is  used  which  is  initially  set  equal  to  zero  on  the  same  ex¬ 
tended  lonum.  The  following  program  computes  one  iteration  of  the  geometric  filter. 


1 . 

a  -- 

l  ,  3  * 

-  0  ,  c  - —  3  ,  d 

—  1. 

2  . 

9  l  71  , 

n »  — 

ma  x  1  f  ( m ,  n ) 

,  mm 

[ f (m-a ,  n-b) -1 , 

f  (m, 

n)  +  11 

l- 

3 . 

t  .  T.  , 

n )  —  max  g  (m,  n)  , 

min 

[g  (m  -  a,  n  -  b) 

,g  (m, 

n )  +  1 

t 

f  j  r 

l  <  m 

<  M,  1  <  n  < 

N. 

4  . 

I  c  i 

“  i  ; 

a  - - a  ,  b  — 

-b,  d 

—  0 ,  go  to  2  . 

If  d 

=  0;  d  ■ 

— 

5  . 

3  '  71  , 

h)  — 

m  l  n  j  f  ( m ,  n ) 

,  max 

[  f  ( m  -  a ,  n  -  o )  +  l 

,  f  (m, 

n) 

tor 

I  <  m 

<  M,  1  <  n  S 

N. 

6  . 

f  i  m , 

hi  •— 

m  n  j  g  ( m ,  n )  , 

max 

(g  (m  -  a ,  n  -  b) 

,  g  (m 

,  n)  - 

i, 

fur 

L  <  m 

S  M ,  1  <  n  < 

N. 

7  . 

It  .1 

=  1: 

a  —  -a ,  b  — 

-b,  d 

—  0 ,  go  to  5 . 

If  a 

=  0, 

d  1 ,  go  to 

8  . 

8. 

If  c 

=  3; 

a  —  0  ,  b  — -  1 

,  c  •— 

2,  go  to  2. 

If  c 

=  2; 

a  —  1 ,  b  —  1 

,  c  — 

1,  go  to  2. 

I  t  c 

=  1; 

a  —  1 ,  b  — -  - 

1 ,  c  —  0 ,  go  to  2  . 

If  c 

=  0; 

stop. 

Example 


Figure  1  is  a  3AR  image  of  an  airport  in  Windsor,  Ontario.  It  was  made  by  the  STAR-1 
airoorne  SAR  system  designed  and  developed  by  the  Environmental  Research  Institute  of 
Michigan  lEKiM).  The  resolution  is  6  m  x  6  m  and  the  pixel  spacing  is  8.4  m.  Figure  2 
snows  tne  result  of  applying  5  iterations  of  the  geometric  filter. 


Note  tnat  small  strong  returns  retain  their  sharp  edges  in  the  filtered  image.  Also 
medium  contrast  edges  retain  their  sharpness  and  low  contrast  edges  are  still  visible  in 
tne  filtered  image.  Figure  3  shows  3-D  plots  of  the  boxed  small  strong  return  in  the  orig¬ 
inal  and  filtered  images.  Figure  4  shows  similar  3-D  plots  for  the  boxed  medium  contrast 
edge . 
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In  ordec  to  measure  tne  "amount"  of  speckle  reduction,  we  define  a  speckle  index  as 
follows.  In  view  of  the  multiplicative  nature  of  speckle  noise  (Goodman^ ,  2 5  •  ,  the 

ratio  of  its  deviation  to  its  mean  seems  to  be  a  reasonaole  measure  of  the  amount  if  5n»r- 
kle  noise  present.  For  1  <  m  £  M  and  1  <  n  £  N,  we  define  an  approximation  to  the  local 
deviation  Dy 

o(m,  n)  =  max  f(m+a,  n+b)-  min  f (m  +  a ,  n  +  b)  .  (1) 

-1  <  a,  d  S  1  -1  <  a,  b  <  1 

The  local  mean  is  defined  as 

1 

„  ( m ,  n )  =  ~  ^  f  ( m  a ,  n  +  b )  <2\ 

a ,  b* - 1 

The  specxle  index  is  then  defined  by, 

^  MN  f 

m=l  n*l 


The  speckle  inaex  for  tne  original  image  shown  in  Figure  1  is  1.05.  The  filtered  male 
shown  in  Figure  2  has  a  speckle  index  of  0.36. 

Synthetic  imagery 

In  order  to  create  an  image  with  synthetic  speckle,  we  begin  by  choosing  a  real-valued 
image,  r(m,  n),  containing  some  patterns  of  interest.  A  random  phase  image,  :in,  n),  is 
generated  where,  for  each  point  (m,  n),  ;(m,  n)  is  an  independent  sample  from  the  uniform 
distribution  over  the  interval  from  0  to  2-r.  We  then  define  a  synthetic  complex  reflectiv¬ 
ity  function  by 

g(m,  n)  =  r (m,  n)  exp  []$(m,  n) ]  <4i 

where  j  *  (-1)1/2.  An  impulse  response  function  is  defined  by 

i(sinc  j  m)  (sine  n)  for  1  m  !  ,  ]n!  S  2, 

i  z  ) 

0  otherwise. 

The  synthetic  complex  radar  image  is  defined  by  h  =  k  *  g  where  *  denotes  convolution. 
Finally,  the  synthetic  detected  image  is  defined  by  f(m,  n)  «  h(m,  n)  .  This  d»tect°d 
image  is  then  modified  by  a  two  step  process  which  is  used  at  ERIM  on  real  radar  imao^rv 
to  reduce  16-bit  signed  data  to  8-bit  unsigned  data  for  display  on  a  CRT.  It  consists  of 
a  linear  scaling  followed  by  a  non-linear  mapping. 

The  image  r(m,  n)  used  to  create  the  synthetic  image  in  Figure  5a  has  a  background  level 

of  12.  It  has  five  rows  of  squares  of  sizes  32  x  32,  16  x  16,  8  x  8,  4  x  4  and  2x2  pix¬ 
els.  The  pixel  values  of  the  squares,  going  from  left  to  right,  are  170,  130,  90  and  50. 

The  two  large  gratings  consist  of  bright  and  dark  stripes  of  dimension  64  x  h  pixels.  The 

bright  stripes  have  a  pixel  value  of  170  in  the  first  grating  and  90  in  the  second.  The 

dark  stripes  are  at  tne  background  level  of  12.  The  stripes  in  the  four  smaller  qratnn'= 

have  dimension  32  x  4.  The  bright  stripes  have  pixel  values  170,  130,  90  and  50,  and  the 

dark  stripes  are  at  the  background  level  of  12. 

Comparison  with  look-averaging 

Look  averaging,  the  noncoherent  addition  of  multiple  images,  is  a  commonly  used  tech¬ 
nique  to  reduce  speckle  noise  in  SAR  imagery  and  is  an  obvious  reference  for  comparina  the 
effectiveness  of  the  geometric  filter.  In  order  to  make  a  fair  comparison  between  the  geo¬ 
metric  filter  and  look-averaging,  the  same  information  was  used  for  both  methods. 

Look-averaging  was  carried  out  as  follows.  The  Fourier  transform  of  the  complex  imaae 
(  *  h(m,  n)  -  see  the  preceeding  section)  was  taken  and  its  square  domain  was  divided  into 
four  smaller  squares.  Each  of  these  four  parts  of  the  Fourier  transform  were  then  inverse- 
transformed  to  obtain  four  complex  looks.  The  detected  looks  were  obtained  by  computing 
the  magnitude  of  the  complex  looks.  Finally,  the  average  of  these  four  detected  looks  was 
computed . 
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Tne  resale  of  tins  look-averaging  process  is  the  image  in  Figure  5b.  The  result  of  f : 
iterations  of  tne  geometric  filter  applied  to  the  original  image  is  in  Figure  5c.  The 
speckle  indices  for  these  tnree  images  are  given  in  Table  1. 

Tab le  1.  Speckle  Indices 


Speckle  Index 


Or  lginal 
Look-average 
Geometric  filter 


0.524 

0.223 

0.065 


Figure  6  snows  3-D  plots  of  the  boxed  square.  Figure  7  shows  2-D  plots  of  a  slice 
through  tne  same  square  and  Figure  8  shows  2-D  plots  of  a  slice  across  part  of  one  of  th« 
gratings. 

It  appears  from  the  above  that,  at  least  in  this  case,  the  geometric  filter  outperforms 
look-averaging.  Thus,  it  could  be  used  either  to  produce  higher  quality  imaaery  or  portae 
to  proouce  imagery  of  tne  same  quality  at  a  lower  cost  from  less  data. 

Summary 

The  geometric  filter  has  become  a  standard  tool  at  ERIM  for  use  in  processmq  SAR 
imagery.  It  iias  oeen  found  to  be  useful  in  preparing  imagery  both  for  human  inspect’."'" 
and  for  computer  algorithms  such  as  edge  detection.  It  reduces  speckle  effectively  while 
at  tne  same  time  preserving  the  spatial  information  in  the  image  such  as  edq°s,  stronq 
returns,  etc.  This  filter  appears  to  outperform  the  commonly  used  look -aver  a g  i  ng  -ethri 
tor  reducing  speckle. 
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Geometric  filter  for  reducing  speckle 


Thomas  R.  Crimmins  Abstract.  A  nonlinear  speckle  filter  based  on  geometric  concepts  s  ce- 

Environmental  Research  Institute  fined,  and  an  example  of  its  effectiveness  on  synthetic  aperture  raaar 

of  Michigan  imagery  is  shown.  Comparison  with  look-averaging  is  made  using  arti- 

Radar  Division  ficial  imagery  with  synthetic  speckle. 
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1.  INTRODUCTION 

The  presence  of  speckle  in  imagery  produced  with  coherent 
illumination  reduces  the  detectability  of  objects  m  the 
images. It  also  reduces  the  effectiveness  of  some  computer 
algorithms  ie.g..  edge  detection)  designed  for  automatic  im¬ 
age  analysis. h 

The  geometric  filter  was  designed  to  reduce  speckle  in 
synthetic  aperture  radar  (SAR)  imagery  while  preserving  the 
spatial  information  in  the  image  such  as  edges,  strong  returns, 
etc.  It  is  an  iterative  nonlinear  algorithm. 

2.  DEFINITION  OF  THE  GEOMETRIC  FILTER 

The  geometric  filter  algorithm  is  based  on  applying  a  single 
iteration  of  an  iterative  convex  hulling  algorithm  alternately 
to  the  image  and  to  its  complement  (negative  of  the  imagei.  It 
is  essentially  a  one-dimensional  algorithm  that  is  applied  suc¬ 
cessively  in  four  different  directions  in  the  two-dimensional 
image:  horizontal,  vertical,  and  the  two  diagonal  directions. 

We  will  describe  the  application  of  the  one-dimensional 
algorithm  in  the  horizontal  direction.  The  image  is  assumed  to 
have  integer  grav-level  values  between  0  and  255. 

Each  horizontal  row  of  pixels  defines  a  one-dimensional 
function  The  "discrete  graph"  of  this  function  can  be  repre¬ 
sented  by  a  subset  of  a  discrete  grid  256  points  high  and  as 
many  points  wide  as  there  are  pixels  in  the  horizontal  direc¬ 
tion  of  the  image  The  umbra  of  the  function  consists  of  all 
points  in  this  grid  on  or  below  its  discrete  graph.  One  iteration 
of  a  convex  hulling  algorithm  is  applied  to  this  umbra  (If 
many  iterations  were  performed,  an  approximation  to  the 
convex  hull  of  the  umbra  would  be  formed  i  Next,  one  itera¬ 
tion  of  this  convex  hulling  algorithm  is  applied  to  the  comple- 

Inv ited  Paper  SO-108  received  Oct  30.  1085,  revised  manuscript  received 
Ian  ".  1986.  accepted  for  publication  Jan  ",  1986.  received  by  Managing 
Editor  Feb  3.  1986  This  paper  is  a  revision  of  Paper  556-30  which  was 
presented  at  the  SP1E  Internationa)  Conference  on  Speckle.  .Aug  30-33. 
1985  San  Diego.  Calif  The  paper  presented  there  appears  lunrefereedi  in 
SPIE  Proceedings  Vol  556 
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ment  of  the  set  resulting  from  the  prev  ious  step  The  output  i' 
the  function  whose  umbra  is  the  complement  ot  the  set  result¬ 
ing  from  this  last  step 

The  procedure  described  above  is  now  applied  in  each  or 
the  other  three  directions  This  constitutes  one  iteration  of  the 
geometric  filter 

For  a  more  detailed  explanation  ot  the  geometric  derr.  a- 
tion.  see  Ret  ”  by  Crimmins  A  program  that  computes  one 
iteration  of  the  geometric  niter  is  presented  ,n  Sec  s  'the 
Appendix)  This  program  was  implemented  on  a  De.An/u  IP 
5500  digital  video  processor  using  a  \A,\  I  I  "so  as  j  host 
system.  Running  time  is  about  14  s  per  iteration  tvir  5  I  2  -512 
images  Usually  from  4  to  10  iterations  are  used 

3.  EXAMPLE 

Figure  1  is  a  SAR  image  of  an  airport  in  Windsor.  Ontario  It 
is  a  four  look  average*  made  by  the  STAR- 1  airborne  SAR 
system  designed  and  developed  by  the  Environmental  Re¬ 
search  Institute  of  Michigan  I  ERIM)  The  resolution  is  b  m  • 
6  m.  and  the  pixel  spacing  is  8.4  m 

Figure  2  shows  the  result  of  applying  live  iterations  of  the 
geometric  filter.  Note  that  small  strong  returns  retain  then- 
sharp  edges  in  the  filtered  image.  Also,  medium-contrast 
edges  retain  their  sharpness,  and  low -contrast  edges  are  'till 
visible  in  the  filtered  image.  Figure  3  'hows  3-D  plots  ot  the 
boxed  small  strong  return  in  the  original  and  filtered  images 
Figure  4  shows  similar  3-D  plots  for  the  boxed  medium- 
contrast  edge 

For  measuring  the  "amount"  of  speckle  reduction,  the 
speckle  index  is  defined  in  the  following  way:  In  view  of  the 
multiplicative  nature  of  speckle  noise  (Goodman. g  p  25 1.  the 
ratio  of  its  deviation  to  its  mean  seems  to  be  a  reasonable 
measure  of  the  amount  of  speckle  noise  present  For  I  s  m  s 
M  and  1  <  n  £  N.  we  define  an  approximation  to  the  local 
deviation  by 

in m.  ni  =  max  tim  -  a.  n  -  bi 

-  I  bs  I 

min  fim  -  a.  n  -  hi  .  ! 

- !  --J  bs  i 

where  f  is  the  function  representing  the  image  The  local 
mean  is  defined  as 

M-im.  ni  =  —  ^  Km  -  a.  n  -  hi  3' 

9  w 

J  6=  -  I 

The  speckle  index  is  then  defined  by 
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Fig.  1.  Original  SAR  image. 
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Fig.  2.  Result  of  five  iterations  of  geometric  filter. 


i  M  * 

I  ^  ^  trim,  ni 

speckle  index  ■  -  >  ?  - 

MN  ^  ^  n.m.  m 

Tl  =  H  =  ■ 


The  speckle  index  tor  the  original  image  shown  in  Fig  I  is 
1  05  The  tillered  image  shown  in  Fig.  2  has  a  speckle  index 
of  0  36. 


4.  SYNTHETIC  IMAGERY 
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Fig.  3.  Three-dimensional  plots  of  a  small  strong  return 
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Fig.  4.  Three-dimensional  plots  of  a  medium-contrast  edge 


interest  A  random  phase  image  6im.  ni  is  generated  where, 
for  each  point  im.  n).  d>im.  n)  is  an  independent  sample  irom 
the  uniform  distribution  over  the  interval  trom  "  to  2-  x\  e 
then  derine  a  synthetic  complex  reflectivity  function  ny 


To  create  an  image  with  synthetic  speckle,  we  begin  by  chocs-  kim.  ni  = 
ing  a  real-valued  image  rtm.  n)  containing  some  patterns  of 


gim.  ni  =  rim.  ni  exp(j<tnm.  ni|  .  a 

where  j  =  i  -  I )'  An  impulse  response  function  is  defined  by 

kim.  ni  =  /  (  smc  ~T~m  )  (  sine  ~n  )  for  m.  n  <  2  . 


0  otherwise 


'^•V  v ■>"'/*  ,  \‘2-“ '■aiftfo,  r.’A'-.A-’.  ■' 


GEOMETRIC  FILTER  FOR  REDUCING  SPECKLE 


(c) 


Fig.  S.  (a)  Original  imaga  with  synthetic  speckle,  (b)  Average  of  four 
looks,  (c)  Result  of  five  iterations  of  geometric  filter. 


The  synthetic  complex  radar  image  is  defined  by  h  =  k  *  g. 
where  *  denotes  convolution.  Finally,  the  synthetic  detected 
image  is  defined  by  f(m.  n)  =  |h(m.  n)|. 

The  image  rim.  n)  used  to  create  the  synthetic  image  in  Fig. 
5(a)  has  a  background  level  of  12.  It  has  five  rows  of  squares 
of  sizes  32*32.  16*16.  8*8.  4*4.  and  2*2  pixels.  The 
pixel  values  of  the  squares,  going  from  left  to  right,  are  170. 
130.  90.  and  50.  The  two  large  gratings  consist  of  bright  and 
dark  stripes  of  dimension  64  <  6  pixels.  The  bright  stripes 
have  a  pixel  value  of  170  in  the  first  grating  and  90  in  the 
second.  The  dark  stripes  are  at  the  background  level  of  12. 
The  stripes  in  the  four  smaller  gratings  have  dimension  32*4. 
The  bright  stripes  have  pixel  values  170.  130.  90.  and  50.  and 
the  dark  stripes  are  at  the  background  level  of  12. 

5.  COMPARISON  WITH  LOOK-AVERAGING 

Look-averaging,  the  noncoherent  addition  of  multiple  statis¬ 
tically  independent  images,  is  a  commonly  used  technique  to 
reduce  speckle  noise  in  SAR  imagery  and  is  an  obvious  refer¬ 
ence  for  comparing  the  effectiveness  of  the  geometric  filter. 
To  make  a  fair  comparison  between  the  geometric  filter  and 
look-averaging,  the  same  information  was  used  for  both 
methods. 

Look-averaging  was  carried  out  as  follows.  The  Fourier 
transform  of  the  complex  image  [h(m.  n) — see  the  preceding 
section]  was  taken,  and  its  square  domain  was  divided  into 
four  smaller  squares.  Each  of  these  four  parts  of  the  Fourier 
transform  were  then  inverse-transformed  to  obtain  four  com¬ 
plex  looks  The  detected  looks  were  obtained  by  computing 
the  magnitude  of  the  complex  looks.  Finally,  the  average  of 
these  four  detected  looks  was  computed. 

The  result  of  this  look-averaging  process  is  the  image  in 
Fig.  5(b).  Figure  5(c)  is  the  result  of  five  iterations  of  the 
geometric  filter  applied  to  the  original  image  The  speckle 
indices  for  these  three  images  are  given  in  Table  I. 

Figure  6  shows  3-D  plots  of  the  boxed  square  selected  from 
Figs.  5(a)  through  5(c).  shown  top  to  bottom,  respectively. 


TABLE  I.  Speckle  Indices  (or  the  Image*  of  Fig.  5 


Image 

Speckle  rnaex 

Original 

0.524 

Look-average 

0.223 

Geometric  filter 

0.065 

Figure  7  shows  2-D  plots  of  a  slice  through  the  same  square, 
and  Fig.  8  shows  2-D  plots  of  a  slice  across  part  of  one  of  the 
gratings. 

It  appears  from  the  above  that,  at  least  in  this  case,  the 
geometric  filter  outperforms  look-averaging.  Thus,  it  could  be 
used  either  to  produce  higher  quality  imagery  or  perhaps  to 
produce  imagery  of  the  same  quality  at  a  lower  cost  from 
fewer  data. 

6.  SUMMARY 

The  geometric  filter  has  become  a  standard  tool  at  ER1M  for 
use  in  processing  SAR  imagery  .  It  has  been  found  to  be  uselul 
in  preparing  imagery  both  for  human  inspection  and  for  com¬ 
puter  algorithms  such  as  edge  detection  It  reduces  speckle 
effectively  while  at  the  same  time  preserving  the  spatial  infor¬ 
mation  in  the  image,  e.g. ,  edges  and  strong  returns.  This  filter 
appears  to  outperform  the  commonly  used  look-averaging 
method  for  reducing  speckle. 
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8.  APPENDIX 

The  image  is  a  function  f(m.  n)  where  I  £  m  £  \1.  I  <nsN. 
and  the  values  of  f  are  integers  between  0  and  255  A  border 
of  zeros  is  added  to  the  image  so  that  f(m.  nl  is  defined  lor 
OsmsM  -  I  andO  snsN  *  I.  An  auxiliary  image  gim.  m 
is  used  that  is  initially  set  equal  to  zero  on  the  same  extended 
domain.  The  following  program  computes  one  iteration  of  the 
geometric  filter. 

1.  a*—  1 .  b  «—  0.  c  «—  3.  d  «—  I. 

2.  g(m.  n)  «—  max{f(m.  n).  min|f(m-a.  n-bi  -  I .  fim.  ni 
+  I  ] } .  for  1  £  m  s  M .  I  £  n  £  N . 

3.  f(m.  n) «—  max{g(m,  n).  min|g(m  -  a.  n  -b).  g(m.  nl  -  I . 
g(m  +  a.  n  *  b)  +  1]}.  for  I  £  m  £  M.  I  £  n  £  V 

4.  If  d  =  I ;  a  ♦ - a,  b  « - b.  d  ■<—  0.  go  to  2. 

If  d  =  0;  d  «—  I .  go  to  5. 

5.  g(m.  n) «—  min{f(m.  n).  max(f(m  -  a.  n  -  b)  ~  I.  f(m.  ni 
-  1  ]}.  for  I  £  m  £  M.  I  £  n  £  N. 

6.  f(m.  n) «—  min{g(m.  n).  max|g(m  -  a,  n  -  b).  g(m.  ni  - 
I .  gun  +  a.  n  +  b)  -  1  ]}.  for  I  £  m  £  M.  I  £  n  £  V 

7.  If  d  —  1 ;  a  ♦ - a.  b  < - b,  d  «—  0,  go  to  5 

If  d  =  0;  d  «—  I .  go  to  8. 

8.  If  c  =  3;  a  «—  0,  b  «—  I .  c  «—  2.  go  to  2. 

If  c  =  2:  a  «—  I .  b  *—  I .  c  «—  1 .  go  to  2. 

If  c  =  I ;  a  «—  I .  b  ♦ - I .  c  «—  0.  go  to  2. 

If  c  =  0;  stop. 
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Fig.  6.  Three-dimensional  plots  of  boxed  square. 


Fig.  8.  Two-dimensional  plots  of  a  slice  across  part  of  a  grating. 


Fig.  7.  Two-dimensional  plots  of  a  slice  through  a  square. 
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The  image  is  a  function  f(j,  k)  where  1  <  j  £  J.  1  £  k  <  K  and 
the  values  of  f  are  non-negative  integers.  A  border  of  zeroes  is 
added  to  the  image  so  that  f(j,  k)  is  defined  for  0  £  j  <  J  +  1  and 
0  <  k  <  K  +  1.  An  auxiliary  image  g(j,  k)  is  used  that  is  initially 
set  equal  to  zero  on  the  same  extended  domain. 


The  control  parameter  n  is  a  positive  integer  and  is  an  input 
parameter.  The  following  program  computes  one  iteration  of  the  gen¬ 
eralized  geometric  filter. 

1.  a  ♦  1,  b  «  0,  c  «  3,  d  «  1. 


2.  g(j,  k)  «  maxjf ( j ,  k),  min[f(j  -  a,  k  -  b)  -  n, 

f(j,  k)  +  n][ ,  for  1  £  j  £  J,  1  £  k  £  K. 

3.  f(j,  k)  <e  max{g(j,  k),  min[g(j  -  a,  k  -  b), 

g(j,  k)  +  n,  g(j  +  a,  k  +  b)  +  n]}, 
for  1  <  J  <  0.  1  £  k  £  K. 

4.  If  d  «  1;  a  «  -a,  b  ♦  -b,  d  «  0,  go  to  2. 

If  d  »  0;  d  «  1,  go  to  5. 

5.  g(j,  k)  <■  min {f ( j ,k ) ,  max[f(j  -  a,  k  -  b)  +  n, 

f ( j »  k)  -  n]f,  for  1  £  j  £  J,  1  £  k  £  K. 


6.  f ( j  ,k)  «  min jg( j ,  k),  max[g(j  -  a,  k  -  b), 

g(j,  k)  -  n,  g( j  +  a,  k  ♦  b)  -  n]}, 
for  1  £  j  £  J,  1  <  k  <  K. 


7.  If  d 
If  d 

8.  If  c 
If  c 
If  c 
If  c 


I; 

a  «  -a 

,  b  «  -b,  d  «  0,  go  to 

0; 

Q. 

* 

M 

go  to  8. 

3; 

a  ♦  0, 

b  ♦  1,  c  *  2,  go  to  2 

2; 

a  «  1, 

b  «  1,  c  *  1,  go  to  2 

1; 

a  «  1, 

b  «  -1,  c  «  0,  go  to  ; 

0; 

stop. 

